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When double-stranded DNA molecules are heated, or exposed to denaturing agents, the two 
strands get separated. The statistical physics of this process has a long history, and is commonly 
described in term of the Poland-Scheraga (PS) model. Crucial to this model is the configurational 
entropy for a melted region (compared to the entropy of an intact region of the same size), quantified 
by the loop factor. In this study we investigate how confinement affects the DNA melting transition, 
by using the loop factor for an ideal Gaussian chain. By subsequent numerical solutions of the 
PS model, we demonstrate that the melting temperature depends on the persistence lengths of 
single-stranded and double-stranded DNA. For realistic values of the persistence lengths the melting 
temperature is predicted to decrease with decreasing channel diameter. We also demonstrate that 
confinement broadens the melting transition. These general findings hold for the three scenarios 
investigated: namely 1. homo-DNA, i.e. identical basepairs along the DNA molecule; 2. random 
sequence DNA, and 3. “real” DNA, here T4 phage DNA. We show that cases 2 and 3 in general 
give rise to broader transitions than case 1. Case 3 exhibits a similar phase transition as case 2 
provided the random sequence DNA has the same ratio of AT to GC basepairs. A simple analytical 
estimate for the shift in melting temperature is provided as a function of nanochannel diameter. 
For homo-DNA, we also present an analytical prediction of the melting probability as a function of 
temperature. 


I. INTRODUCTION 

The three dimensional structure of double-stranded 
Deoxyribonucleic acid (dsDNA) DNA is the famed 
Watson-Crick double-helix, within a broad range of salt 
and temperature conditions.® This double-helix encodes 
the genetic information of an entire organism in terms 
of four types of bases. Base complementarity guarantees 
that the same information is contained in the two strands 
of a dsDNA.® 

By temperature increase the double-stranded double¬ 
helical DNA progressively denatures - DNA melting, see 
Figure Partially melted DNA is an alternating set of 
intact, stiff, double-stranded regions (persistence length 
about 50 nm) and of single-stranded, floppy (persistence 
length about 1 to a few nm), “melted” regions (DNA 
bubbles).® Each of these regions is subject to thermal 
fluctuations. Interestingly, the average free energy as¬ 
sociated with breaking an AT-basepair is smaller than 


the corresponding energy for a GC-basepair.® In melting 
studies, therefore, DNA tends to melt first in AT-rich 
regions (45.79 °C at salt concentration 0.01 M for pure 
AT) and only at higher temperature (95.88 °C at salt 
concentration 0.01 M for pure GC regions) in GC-rich 
regions. 

On the theoretical side. The Poland-Scheraga (PS) 
model of DNA melting has been proven to well repro¬ 
duce (macroscopic) melting data.®^®2l Note, however, 
that there is an alternative model, the Peyrard-Bishop 
model,li®i^ see Ref. |7| for a comparative study of the 
two models. Herein we use the PS model, which is 
an Ising model with a long-range term due to the en¬ 
tropy associated with the melted single-stranded regions, 
and has the following parameters: two hydrogen bond 
(Watson-Crick) energies (AT or GC bonds), ten inde¬ 
pendent stacking (nearest neighbor) parameters, the loop 
factor (or, loop functioiP), and the ring-factor (bub¬ 
ble initiation) parameter (related to the cooperativity 
parameter).® Rather recently the ten stacking parame- 
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FIG. 1. Cartoon of DNA melting in a nanochannel. A given 
microstate of DNA, at elevated temperatures, is an alternat¬ 
ing set of intact double-stranded regions and melted single¬ 
strand regions (DNA bubbles). Within the Poland-Scheraga 
model, the statistical physics of DNA melting is described in 
terms of two hydrogen-bond energies, ten stacking (nearest 
neighbor) parameters and the loop factor g(rn) which quan¬ 
tifies the conformational entropy for a melted region corre¬ 
sponding to m basepairs (relative to the entropy of an un¬ 
melted region of the same size). The functional form of g{m) 
is different for unconfined and nanoconfined DNA. 


ters were independently measured for the first time for 
different temperatures and salt concentrations.^ 

Theoretical interest in the melting transition origi¬ 
nates, in part, from the fact that the PS model exhib its a 
phase-transition which can be analyzed exactly.^SEH This 
analysis treats unconfined DNA molecules, with identical 
hydrogen bond and stacking energies (homo-DNA). For 
unconfined DNA the loop factor g{m) (see above) scales 
with DNA bubble size to as a power-law: g(m) ~ m~‘^ 
with a loop exponent c. The exponent c determines the 
order of the phase transition; if c > 2 the melting transi¬ 
tion is first order 

More recent studieJ^^^^ address the challenge of un¬ 
derstanding how the DNA sequence affects the melting 
transition, typically by analyzing the melting of random 
sequence DNA.122H211 Unlike for the case of homo-DNA 
melting, no exact results exist for these kinds of systems. 

Herein, we go beyond previous studies by considering 
how confinement affects the DNA melting transition. In 
particular, we seek to quantify how confinement affects 
the DNA melting temperature and the width of the melt¬ 
ing transition. In contrast to Ref. [571 where molecular 
dynamics simulations of very short DNA (15 basepairs) 
were performed, we are here interested in the thermo¬ 
dynamic limit (we study typical DNA sizes of 2 • 10® 
basepairs). In Ref. [551 we calculate how confinement in¬ 
fluences 5 (to) and show how this influences melting for 
a highly simplified model system. In the present study, 
we go beyond Ref. [55] by calculating what consequences 
the effect of confinement on g{m) has for the melting 
of a more realistic model of DNA. Using as input the 
functional form for g{m)^^ we provide results of exten¬ 
sive simulations within the Poland-Scheraga model for 
I. homo-DNA, 2. random sequence DNA and 3. “real” 
DNA (T4 phage). We also provide a simple formula for 
estimating how the melting temperature changes as a 
function of the channel diameter and present an analyt¬ 
ical prediction of the melting probability of homo-DNA 


as a function of temperature. 

The present study is inspired by recent experiments 
which demonstrate the potential to study local proper¬ 
ties of the DNA melting transition in nanochannels.^^SMl 
Other recent experime nts no w allow studying also the 
DNA melting dynamics.^^^^ 

The biological relevance of DNA melting in confined 
environments is, for instance, due to the fact that local 
denaturation of DN A is n ecessary for protein binding to 
a DNA single- stran d,and is implicated in transcrip¬ 
tion initiation.lS12SI 


II. REVIEW OF THE POLAND-SCHERAGA 
MODEL 

In this section we introduce the Poland-Scheraga (PS) 
model for unconfined DNA. In the next section we discuss 
how the PS model must be modified in order to study 
DNA melting under confinement. 


A. General considerations 


Consider double-stranded DNA with N internal base- 
pairs, that are clamped at both ends for simplicity 
(Fig.i. The free energy required for breaking a hydro¬ 
gen bond at basepair i is denoted by i?hb(*)- There are 
two values for this parameter (AT-bonds or GC-bonds). 
For disrupting the stacking (nearest neighbor) interac¬ 
tions between basepairs i — 1 and i there is a free energy 
cost Ast(* — 1, i). There are ten stacking parameters.!^ In 
addition we need the ring factor ^ « I0“® (see Ref. dsi) 
which is a Boltzmann factor associated with the free en¬ 
ergy cost of creating two “boundaries” between intact 
and melted DNA. Note that E\±,{i) and Est{i — 1, i) have 
energetic as well as entropic contributions: £'hb(*) = 
t^hb(*)-'C5'hb and E^tii-lfi) = Ustii-lfi)-TSst, where 
U denote energetic contributions, T is the temperature 
and S'hb and Sst are entropies, commonly assumed to 
be basepair independent. Finally, the Poland-Scheraga 
model requires as input the ratio of the number of con¬ 
figurations for a melted region (DNA bubble) and the 
number of configurations of an intact region, quantified 
by the loop factor g{m), see sections III A and IIIB 


B. The melting temperature and the 
Marmur-Doty formula 

Throughout the main text we repeatedly refer to the 
melting temperature, Tm, or deviations therefrom. Tm is 
defined as the temperature where the melting probability 
= 1 / 2 . 

We can estimate the melting temperature, Tm by set¬ 
ting the total free energy per basepair, E\^\^ -\- Est, equal 
to zero. This leads to 
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Tm = (C/hb + U,t)/Srei (1) 

with Srei = 5'hb + *S'st, and we find that the free energy 
can also be expressed as 


E = S,ei{.TM-T). ( 2 ) 

The total free energy of a given DNA sequence can be 
estimated as the sum of the free energies of the AT and 
the GC portion: 

E = /atAAat + /gcAAgc 

= <5'ref(/AT'rAT + /gC?GC - T) (3) 

with the free energy of random AT and random GC se¬ 
quences AAat and AEqc- 

With Eq. ([^ we find the Marmur-Doty formula, which 
gives a simple estimate for the melting temperature, Tm 
of unconfined DNA: 


Tmd = /at "Tat + focToc 


(4) 


where /at/gc is the fraction of AT/GC bonds in the 
DNA sequence. Tat/gc is the melting temperature for 
pure random AT/GC sequences. We note that this for¬ 
mula is rather crude, in reality Tm depends weakly also 
on ^ and g{m). 

We will later use the stability parameters introduced 
in Sec. IIA from Ref. [T3I (Table 1). Herein, the ten 


stacking parameters were experimentally determined and 
the hydrogen bond energies were found from the empiri¬ 
cal Marmur-Doty melting temperature together with Eq. 
([^ and 


= AEat - 4 


E 




(5) 


at,ta,aa,tt 


and the corresponding equation for GC. Note that we 
use a different notation than as Ref. [131 E = —AG and 
Eref = -AS. 

Eor homo-DNA, the melting temperature is estimated 
as: 


rrAA , r/AT 
nA ^ + '^hb 

^ ^hb + S',* 

= Tat + 


( 6 ) 

(7) 


At salt concentration O.OI M, we have = 52.73°C. 

Even though investigation of unconfined DNA is not 
the main purpose of this study, we elaborate a bit on de¬ 
viations from the Marmur-Doty formula in appendix |C| 
where we also investigate the properties of the DNA melt¬ 
ing temperature on the value of c for unconfined DNA. 


III. THE LOOP FACTOR 


A. Loop factor for unconfined DNA 

Let us now consider g{m) for unconfined DNA (see 
also Ref. 1^51) . For simplicity, assume that the melted 
DNA region is described by a self-avoiding random walk 
of 2m steps in 3 dimensions on a lattice. If we denote by z 
the coordination number (number of nearest neighbors in 
the absence of self-exclusion), and limit ourselves to large 
number of steps, the number of configurations for a closed 
loop (ring polymer) is Dj-ing = where c = 1.76. 

For a linear polymer (intact DNA), similarly Diinear = 
5™, with a different coordination number z, since the per¬ 
sistence length is in general different for double-stranded 
DNA compared to single-stranded DNA. The prefactor 
(z/z)^"* is included in formalism by a redefinition of the 
the Watson-Crick energy Ehb —t Ehb -k 2i?Tlog(z/z^/^) 
and the loop factor then becomes 

gu{m) ~ (8) 

for unconfined DNA. If one includes self-avoidance be¬ 
tween intact DNA regions and melted regions one finds 
c = 2.1^. For an unconfined ideal (phantom) chain one 
obtains c = 3/2. A few words of caution are here needed. 
Strictly speaking, the quantity g{m) is proportional to 
the number of configurations of a random walk of 2m 
steps with the same start and end point, where the ran¬ 
dom walk is constrained to have no bound complementary 
basepairs along its path.^ We here follow the approxima¬ 
tion common to the literature,^ namely, assuming that 
g{m) can be estimated by removing the constraint of no 
bound states. In Ref. HSjwe discuss the effects caused by 
including the constraint mentioned above. 

The power-law form for gii(m) gives rise to an 
effectiv ely lo ng-range interaction for molten DNA 
regions Note that the persistence length of single- 
stranded DNA is 1-5 nm. So, the loop-correction g{m) ^ 
should not be included for small bubbles, although 
it is common practice to include it for all m-values. A 
weight of ^g{m) is assigned to each bubble. Therefore 
the value of the ring factor ^ from Ref. [T31 which was 
obtained from very small DNA molecules neglecting the 
loop-correction, has to be complemented by a loop fac¬ 
tor ( 7 ii(m) = m~‘^. We note that a useful form for g(m) 
approximately correct for all m is given in, for instance. 
Ref. IHl 


B. Loop factor for DNA confined in a nanochannei 

The aim of this study is to understand how chan¬ 
nel confinement influences the melting transition for real 
DNA. To that end we need to find a functional form for 
g{m) for a polymer confined to a nanochannel. In general 
we have: 


In this section we give details about the loop factor for 
confined and unconfined DNA molecules. 


gc{m) 


GssjNss) 

Gds{Nds) 


(9) 
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The quantity G'ds(-^ds) is Green’s function for a double- 
stranded DNA consisting of N^s Kuhn lengths in a square 
channel with side length D and random initial positions 
for the polymer start position. Similarly, Gss(Ass) is pro¬ 
portional to Green’s function for a single-stranded DNA 
region of Agg Kuhn lengths with identical start and end 
positions [below we determine the constant of propor¬ 
tionality using known asymptotic results for g{m)\. The 
number of base pairs is denoted by m and related to 
the Kuhn lengths according to Agg = 2m/{x£ss) and 
Nds = fn/{x£ds) with a conversion factor x = 3bp/nm 
(the center-to-center distance between adjacent basepairs 
is 0.34 nm, see Ref. HI). We use £ss = 6nm for single- 
stranded DNA and ^ds = lOOnm for double-stranded 
DNA throughout this study. Determining Gds(Ads) and 
Gsg(Ags) poses the formidable, unsolved, task of deriv¬ 
ing an expression for the number of conformations for 
a worm-like chain modeP^ with two constraints: (i) no 
polymer conformations can “pass” the channel walls, and 
(ii) the polymer conformations cannot self-intersect. In 
order to provide expressions for Ggg(Asg) we neglect self¬ 
avoidance, i.e., the constraint (ii) above. Notice that 
for typ ical experimental channel sizes, D ~ 50 — 100 
nm EH thus Iss/D <C 1. Single-stranded regions are there¬ 
fore well described by the statistics of an ideal Gaussian 
chain, a nd Ggg (Agg) is given by the solution to a diffusion 
equation.lSSEni The result is given in the Supplementary 
Information of Ref. HU where we find that for an ideal 
Gaussian chain we have 


issue we write 


Gds(A'ds) = i{Nds)Gf^{Nds) (12) 


where 


(Adg) = 


8 1 
7r2 ^ {2k -K 1)2 


exp Adg) 




6D2 


(13) 

(14) 


The quantity G[[‘®(Adg) is the Green’s function for an 
ideal Gaussian chain in weak confinement {D 3> ^ds)-^ 
The factor 7 (Adg) quantifies deviations between the true 
statistics for a double-stranded region and the ideal 
Gaussian case. In order to estimate 7 (Adg) we note that 
in Ref. HT] an interpolation formula for the free energy 
of confinement was given which encompasses both the 
Odijk regime and the diffusive regime of an ideal worm¬ 
like chain. We here choose 7 (Adg) so that the free energy 
of confinement, i?log[Gds(Adg)], equals the one given in 
Eq. (13) of Ref. UTJ in the limit Adg —>■ oo. We then have 

7 (Adg) = exp (2Af (1 - K{£ds/D)) Adg), (15) 

n{z) = {l + 1.672Z -I- 1.287z^)~^^^ . (16) 


Ggg(Agg) = AA-G2 


(-Afc'Ags) 


A 


ss 

k 


6D2 


( 10 ) 

( 11 ) 


where A|® are the eigenvalues of the diffusion operator. 
The constant of proportionality, A, is obtained below. 

For deriving an expression for Gds(Adg) we notice that, 
in contrast to iss/D, the ratio of double-strand Kuhn 
length to channel diameter, £ds/D, is not necessarily a 
small number for realistic channels sizes. To address this 


In the derivation of the above expressions, we have as¬ 
sumed that the entropy of an intact or melted section in 
the interior of a polymer is the same as that of an isolated 
polymer of the same size. We briefly discuss the effect 
of this assumption in the Supplemental Material of Ref. 

m 


Our final task is now to determine the constant A in 
Eq. (101. To that end we require that gc{'m) agrees with 
the unconfined case for a ideal chain, guini) = as 

D ^ oo. In order to analyze the large D limit of gc{m), 
Eqs. (10) and (13) are not suitable. We therefore rather 
use the equivalent resummed expressions 


Ggg(Agg) = AA-G2 


I 6112 

e^Ag, 


/ 6D^n^\ 1 

2^ exp-I -H 


4A^ss 


(17) 


Gr (A^ds) = 


f) / -I 

1 — -j= + 2 ^(—1)" |v^exp {—Bn?^ — V^nerfc (^VBn^ | 


with B — 


3D2 


2 ^dsA^ds 


(18) 


where Eq. 0 was obtained using the Poisson resum¬ 
mation formula, and Eq. (181 is derived in Ref. HU (chap. 


I 

3). With limo^oo Gds(Adg) = 1 we find 

^ _ 27r / 2 /£ ^ ^ 

(~D 


(19) 
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With this result, the known value of the ring factor for 
free DNA ^ = 10“^ can be used. 

The above expressions for g{m) for a channel allow us 
to obtain numerical results for the melting behavior. To 
that end we use the Fixman-Freire approximation with 
stability parameters from Ref. see appendix |A 2| for 
details. In practice only the first term of the sums in Eqs. 

- ( [I^ (fc = 1 and n = 1 ) are sufficient (we denote 
them by Gis\ , gH'^ and G^^^ in the following). We 
use 5 c (w) = GssV^ds^ for small bubbles with size m < 
mo, gdm) = G^sV^ds^ for the intermediate range with 

bubbles of size toq < m < mi and 5 c(m) = GssV^ds^ for 
large bubbles with size m > mi. The transition points 
mo and mi are found numerically as mo = 0.7(I?/4s)^ 
and mi = 17.2{D/£ssd with maximum errors 0.016% and 
0.004%. Also, we “pull out” the hrst term in Eq. (10) 


and in Eq. (13) and redefin e the hydrogen bond energies 
accordingly (see subsection III A), i.e., we write 


/.oAds \ 

5c(m) = /(m)exp ( 20 ) 

Using the explicit form of the eigenvalues of the diffusion 
operator, we have 


9c{m) = f{m) exp 


(. 


\3xD^ 


{£dsK{ids/D) - 24s) m 


( 21 ) 

where we introduced a shifted loop factor /(m). Intro¬ 
ducing f{m) and the associated redefinition of the hydro¬ 


gen bond energies (compare to subsection III A), allows 
us to avoid numerical problems with the FF algorithm 
due to multiplication of exponentially small numbers. 
The function /(m) is plotted in Fig. where we see 
that for small bubbles we have /(m) oc whereas 

for large m we have /(m) oc The full functional 

form smoothly interpolates between these two extremes. 
In Ref. HSIwe show that /(m) has a simple interpretation 
as being proportional to the probability that a confined 
polymer forms a loop. 

In the derivation of the loop factor g{m) we make a 
number of simplifying assumptions. In particular, the 
ideal chain approximation is not a very realistic descrip¬ 
tion of single stranded DNA. However, in Ref. [28|we find 
in simulations of a simple model of DNA that the ef¬ 
fect of confinement on the melting transition is qualita¬ 
tively similar for ideal and self-avoiding chains. For this 
reason, we expect that the 5 (to) used herein provides 
qualitatively correct predictions for the effect of channel 
confinement on the melting transition of DNA. 


IV. RESULTS 

The functional form for the loop factor 5 (to) presented 
in the previous section allows us to study how the melt¬ 
ing transition changes as the channel diameter is var¬ 
ied. We limit ourself to long DNA molecules (thermo- 



FIG. 2. Shifted loop factor /(to) for a polymer confined to a 
square channel of diameter D. The black solid curve shows 
the exact result for an ideal chain contained in Eqs. (§-([^. 
In general we identified three regimes separated by transition 
points Too and toi (see main text). For small and large to we 
have that a power-law behavior, /(to) oc and f{m) oc 

respectively. For illustrative purposes we used a small 
value for the ratio of channel diameter and single-stranded 
DNA persistence length, D/4a = 2. 


dynamic limit) and study melting of three types of pro¬ 
totypical DNA: I. Homo-DNA, i.e. all basepairs are 
identical. For this case we obtain analytical estimates 
alongside the numerical results. 2. Random DNA, i.e., 
the A, T, G or C basepairs are chosen with probabili¬ 
ties, Pa/t and Po/c = 1 ~ Pa/t- 3 - Finally, we con¬ 
sider a “real” DNA sequence, namely that of T4 phage 
DNA. The numerical results are obtained by solving the 
Poland-Scheraga model (Fixman-Freire approximation) 
using the functional form for 5 (to) given above. Details 
are found in appendix [Aj 


A. homo-DNA 


We first consider homo-DNA, i.e. assume that the pa- 

are 


rameters i?hb(f) and Egtii — I, i), see subsection IIA 


independent on i. We first present analytical estimates 
for the melting curve, the melting temperature, and the 
width of the melting curve for small channels, before com¬ 
paring to numerical results. 

In the limit of small channels (or large bubbles), Eq. 
© becomes 


lim 5 c (to) 

m—>-oo 


G: 


( 2 ) 


to(2) 

^ds 


X Aexp 


/ 2 to \ 

64 \a;4s J 

(teIU i^dsK{£ds/D) - 24s) TO^ 

( 22 ) 


Using this approximation, we can derive a number of 
exact results. 
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— — - c= 2.12 

— — - c = 1.76 

- fi ee DNA 

- D/k> = 15 

D/l^ = 10 
• analytic 


FIG. 3. Melting probability as a function of temperature for homo DNA at different channel diameters. Numerical results 
(solid lines) and the analytical prediction (red marks) from Eq. (25l. Notice that decreasing the channel diameter leads to a 
decrease in the melting temperature and to a broadening of the melting curves. The analytical prediction agrees rather well 
with the numerical results already for moderate confinement. Dashed lines show melting probabilities for different values of 
the loop exponent c. We see only a very slight chani 
and comparison to the analytic prediction see figs. [D 

input sequence = ‘AAAAA ... AAAA’, and salt concentration = 0.01 M. Note that 
(unconfined ideal polymer). 


:e of the width of the transition, for further details on the effect of c 
and 12 in the appendix. Parameters: number of basepairs = 2 • 10®, 

free DNA’ here refers to the case c = 1.5 


First, there is a channel induced shift in the melt¬ 
ing temperature: as discussed in the previous sec¬ 
tion, the exponential term in Eq. (211 can be in¬ 
cluded in the hydrogen bond energies, Ehb —>■ Ehb — 
RTi£dsK{ids/D) — 24s) in analogy with how we in¬ 
cluded the coordination number in Sec. IIII Al This leads 
to a shift of the melting temperature according to 


Tm. 


_U _ 

confined — — 


with 




Tm 


(23) 


AS = AS,,-ASds = R 


3x1)2 


(4sk(4s/£>) - 24s) (24) 


and where R is the molar gas constant. Confining an 
ideal (linear) chain to a channel decreases the entropy by 
an amount ASss per link for the single-stranded (melted) 
region and AS'us for the intact double stranded region. 
This entropy of confinement agrees with a scaling argu¬ 
ment given in section 1.1.3 in Ref. For i^sKi^ds/D) > 
2iss we predict that there is an entropy driven decrease 
in melting temperature for small channels. In particular, 
we show in Fig. that for typical experimental channel 
sizes, confinement indeed decreases the melting tempera¬ 
ture of DNA. Moreover, we find that the melting temper¬ 
ature as a function of (-ss/D has a minimum at around 
£ss/D « 0.12. Within our model, this finding follows 
from the fact that a double-stranded region enters the 
Odijk regime (where the entropy of con fineme nt scales 
as {(-ds/DY^^ per Kuhn length [see Eqs. ( 13|l6 l], rather 
than as (4s/41)^ as for the Gaussian chain regime) at 
a larger value of D than a single-stranded region. Note 
that we for the single-stranded region actually do not 
include the Odijk regime in our expression for the asso¬ 
ciated Green’s function. Hence, we cannot realistically 


consider larger values for iss/D than those shown in Fig. 
^ We here note that in the simplified model in Ref. 
25] the persistence length of single-stranded and double- 
stranded regions are the same which rather leads to a 
slight increase of the melting temperature for decreasing 
channel diameters. 

Second, once we have redefined the melting temper¬ 
ature, we can in fact predict the full melting curve. 
To that end we define an effective ring factor as 
= (^V64 )(x4s/ 2)1/2H^ = with K = 

7r®/(48x4s); for 4s= 6 nm we have K = 0.354. Then the 
problem at hand becomes that of melting with a loop 
factor ~ (“one dimensional” random walk) and 

ring factor (all valid for small channels). The change 
in the value of the ring-factor, ^ —>■ is again due to 

entropic confinement effects. We here restrict ourselves 
to D values such that D/lss 4 7, for which we have 
< C- For smaller values of values of T)/4s the diffu¬ 
sion approximation used for the single-stranded regions 
breaks down, see previous section. Using standard ana¬ 
lytical techniques (appendix]^ we find that the expected 
fraction of melted basepairs is given by: 


P = 1 - 


1 

1 -I- (To/3-izoLic_i(zo) 


(25) 


where c = Cid = 1/2 and /3 = exp[AS'(T — 

Tm, condaed)/(RT)] with the “universal” dimensionless 
DNA constant AS/ks ~ 12.51 [since R = 1.987 cal/(mol 
K) and the entropy loss is A^ = 24.85 cal/(mol K), 
see Ref. [T3]. The quantity o-q = oxp[~Flst/(G'F)] is 
the cooperativity parameter, here modified to take into 
account effects due to the channel (th roug h 4ff)- The 
quantity zq is obtained by solving Eq. (B9) numerically; 
whenever no root Zq is found in the interval 0 < Zq < Ij 
one sets P = 1. 
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FIG. 4. Confinement induced shift in melting temperature 
(units of Kelvin) and width of melting curves, for homo-DNA, 
as a function of inverse channel diameter. (Top) Shift of the 
melting temperature, where the theoretical prediction is given 
in Eqs. (231 and \2A\ . (Bottom) Width of the meltin g cu rve, 
where the theoretical result is obtained through Eq. (25 I, see 
main text for details. Input parameters were the same as in 

Fig.m 


Let us now present results of numerical simulations. 
In Fig. we display the melting probability for differ¬ 
ent channel diameters as a function of temperature. The 
numerical results are compared to the analytical predic¬ 
tion of Eq. (251. We also include results for unconfined 
DNA with c = 1.76 and c = 2.12. The homo-DNA melt¬ 
ing curves are generated by using as input a sequence 
AAA...AAA to our Fixman-Freire code. We find that as 
the channel diameter is decreased the melting tempera¬ 
ture increases and the width of the transition increases. 
The melting curve is slightly asymmetric in agreement 
with the prediction in Eq. (251. However, the asymme¬ 
try is less pronounced for melting in the channel com¬ 
pared to melting of unconfined DNA (c = 3/2, c = 1.76 
and c = 2.12). In Fig. we also show the analytical 
prediction from Eq. ( |25[ ) and find good agreement with 
the numerical results. As input for the unconfined DNA 
melting temperature, Tm, appearing in Eq. (231, we use 
the temperature at which P — 1/2 from simulations for 
the unconfined case (c = 3/2). Note that Tm obtained in 
this way deviates by roughly 0.03 °C degrees compared 
to the Marmur-Doty formula. 

Figure shows results for the melting temperature 
(top) and width (bottom) of the melting curve as a func¬ 


tion of channel diameter. These results where extracted 
from melting curves like those presented in Figure the 
melting temperature is the temperature at which we have 
P = 1/2 and the width a of the melting transition is ob¬ 
tained as tr = T(/ = 0.97) — T{f = 0.03). We find 
excellent agreement with the analytical prediction [Eqs. 
(23) and (24)] for the decrease in melting temperature 
with decrease in channel diameter. We get our ana¬ 
lytic prediction for cr by setting P = Pupper = 0.97 and 
P = Piower = 0.03 in Eq. (25) and then solve numeri¬ 
cally (bisection method) in order to obtain temperatures 
Tupper and Tiower respectively. From these solutions we 
then calculate a = Pupper — Piower • Our results show that 
stronger confinement leads to an increase in the width 
of the melting transition. At strong confinement we find 
good agreement between the numerical results and our 
analytical prediction. Since we assume .^ss = 6nm, the 
choice of nanochannel diameters considered in Fig. 
corresponds to 40 nm < D < 90 nm. 


B. Random DNA 


In order to investigate how the sequence affects the 
melting transition, in this section we present numerical 
results for confined random DNA (for such a scenario 
no analytical predictions are available). In all simula¬ 
tions we generated 200 random sequences with a given 
AT-fraction, pat, then used the Fixman-Freire approx¬ 
imation to predict a set of probability profiles. Finally, 
we calculated the average melting probability from this 
set. 

In Fig. [^we display the average melting probability 
for different channel diameters as a function of tempera¬ 
ture for three different pat = {0,0.25,0.5}. The results 
from Fig. [^are further analyzed in Figure|^which shows 
results for the melting temperature (top) and width (bot¬ 
tom) of the melting curve as a function of channel diam¬ 
eter. These results where obtained from melting curves 
as displayed in Fig. in an identical fashion as in the 
previous subsection. 

Figs. and show a number of interesting results. 
The main effect of the heterogeneity (non-zero pat) is to 
further broaden the transition - the transition is broadest 
for the case /at = 1/2, and sharpest for homo-DNA, 
PAT = 1 (pure AT) and pat = 0 (pure GC). Just as for 
the homo-cases (see previous subsection) we find, for a 
fixed value of pati that increasing confinement leads to 
a decreased melting temperature and an increase in the 
width of the transition. Interestingly, the shift in melting 


temperature, Eqs. (23) and (24), obtained for a homo- 
DNA scenario, works rather well also for random DNA, 
see Fig. (top). 

Studying Fig. (Bottom) we see that while confine¬ 
ment broadens the melting transition also for random 
DNA, the effect is relatively small compared to the in¬ 
trinsic width of the transition, which is much larger for 
random DNA compared to the homo-DNA case. This 
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FIG. 5. Average melting probability for random DNA (with different AT fractions) as a function of temperature at different 
channel diameters. For a fixed AT fraction, pat, we find that smaller channel diameters leads to decreased melting temperatures 
and a broader transition (compare to Fig. |^. The main effect of the heterogeneity (non-zero pat) is to further broaden the 
transition compared to melting of homo-DNA - the transition is broadest for the case pat = 1/2. Parameters: number of 
basepairs = 2 • 10®, and salt concentration = 0.01 M. Each curve is an average over 200 random sequences of length 200 kilo 
basepairs with a given AT-fraction, pat- Note that ‘free DNA’ here refers to the case c = 1.5 (unconfined ideal polymer). 
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FIG. 6. Confinement induced shift in melting temperature 
(units of Kelvin) and width of melting curves, for random 
DNA with AT fraction pat , as a function of inverse channel 
diameter. (Top) Shift of the melting temperature. Interest¬ 
ingly, the shift in melting temperature, Eqs. (231 and (241, 
obtained for a homo-DNA scenario, works rather well also 
for random DNA. (Bottom) Width of the melting curve. We 
notice that the width increases with decreasing channel diam¬ 
eter, and that the width is at maximum for pat = 1/2. Input 
parameters were the same as in Fig. 


effect is also seen in Fig. where we find that the dif¬ 
ference between “free DNA” and confined random DNA 
is small. 


C. T4 phage DNA molecules 

We now investigate a “real” DNA sequence - T4 phage, 
and compare to the melting behavior for random DNA. 

In Fig. we display the melting probability for free 
DNA and DNA confined to a channel with different di¬ 
ameters as a function of temperature for T4 phage and 
for random DNA with the same length L = 165643 and 
AT ratio pat = 0.64711. Figures (top, bottom) show 
the shift of the melting temperature and the change of 
the width a for different confinement. 

From Figures and a few things are worth pointing 
out. First, we notice that the melting curves for random 
DNA and T4 phage DNA are very similar. The main 
difference is that the width is somewhat increased for T4 
DNA (however, less than 10 % for all channel diameters 
investigated). Obviously, these findings cannot be gen¬ 
eralized to arbitrary “real” sequences, but shall, rather, 
be used as a motivation for why studying random DNA 
may be relevant also for DNA melting experiments. 


V. SUMMARY AND OUTLOOK 


In this study we investigate how confinement affects 
the DNA melting transition. The statistical physics de¬ 
scription of this process has a long history, and is com¬ 
monly described in terms of the Poland-Scheraga (PS) 
model. The PS model is an Ising type model with hetero¬ 
geneous interaction due to the sequence dependent ener¬ 
gies required for breaking the bonds between the comple- 
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mentary bases. Moreover, including the configurational 
entropy associated with a melted DNA region gives rise 
to an effective long range interaction. It has been found 
previously that this loop entropy factor changes the na¬ 
ture of the DNA melting transition - from “smooth” in 
the absence of this factor, to first order when the loop 
factor is included.^ These results pertain to a scenario 
of identical basepairs along an unconfined DNA. 


FIG. 7. Melting probability for a T4 phage DNA, and random 
DNA, as a function of temperature at different channel diam¬ 
eters. The melting behavior of T4 phage DNA (same pat 
and number of basepairs as the T4 phage DNA sequence) 
is very similar to its associated random sequence melting 
curve. However, the random DNA melting curve is some¬ 
what steeper. Parameters: number of basepairs = 168900, 
and salt concentration = 0.01 M. Each random curve is an 
average over 200 random sequences with a given AT-fraction, 
Pat = 0.647016. The sequence for T4 phage was downloaded 
from the NCBI GenBank (NC_000866.4). 



Inspired by recent experimentswe calculate the 
effect that confinement in a nanochannel has on the melt¬ 
ing transition of DNA. By subsequent numerical solutions 
of the PS model, and simplified arguments, we demon¬ 
strate that the confinement has three effects: the melt¬ 
ing temperature is decreased, the confinement broad¬ 
ens the melting transition, and the shape of the melting 
curve goes from a functional form for unconfined melting 
(loop exponent c = 3/2) to a form consistent with one¬ 
dimensional melting (c = 1/2). We further quantify the 
sequence dependence of the melting transition, and find 
that the transition is broadest when the fraction of AT 
basepairs is 1/2. Relative to the width of the transition 
in the unconfined case, the broadening effect of confine¬ 
ment is smaller for random DNA than for homo-DNA. 
We also investigate “real” DNA, here T4 phage DNA, 
and find that its melting behavior is very similar to the 
melting of a random DNA with the same AT fraction. 
This hnding gives further importance to recent efforts 
for understanding DNA melting of unconfined random 
DNA.E3I2I1 


Throughout this study we use a functional form of the 
loop factor g{m) for a confined chain, which is derived as¬ 
suming ideal chain statistics (no self-exclusion effects in¬ 
cluded), and relaxing the constraint that within a melted 
region no complementary base pair binding must take 
place.l^ While these approximations imply that one can¬ 
not expect the predictions of this study to be exactly 
valid for real DNA, we hnd a qualitatively similar effect 
of confinement upon the melting transition for simula¬ 
tions of an ideal DNA model and of a self-avoiding one 
in Ref. [551 We therefore expect that the results of this 
study furnish a qualitatively correct description of the 
effect of confinement on the melting transition of real 
DNA. 


FIG. 8. Shift in melting temperature and width of melting 
curves T4 phage DNA and associated random DNA (same 
AT fraction), see Fig. as a function of inverse channel 
diameter. (Top) Ghange of the melting temperature (units of 
Kelvin)due to confinement, where the solid curve is Eqs. (231 


and ( |24| . Notice that the shift in melting temperature for T4 
phage DNA is very similar to the melting temperature for a 
random DNA with the same AT fraction. (Bottom) Width of 
the melting curve at different confinement. The width of the 
melting curve for T4 phage DNA is slightly larger than the 
width for a random DNA with the same AT fraction. Input 
parameters were the same as in Fig. 


In a broader perspective, our study sheds light on the 
DNA strand separation in cells. The cell interior is a 
confined and crowded environment, where the crowding 
due to other macromolecules is known to affect vital bi¬ 
ological processes such as the folding of proteins, diffu¬ 
sive transport and reaction rates.l^Our study serves as a 
first step towards an understanding of how confinement 
and crowding affects the DNA melting, with potential 
implications for understanding of the proce sses whereby 
proteins access si ngle s trands of DNA,ISE11 and for tran¬ 
scription initiatiorP^ in living cells. 
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Appendix A: Numerical implementation 

In this appendix we introduce one exact scheme, the 
Poland algorithm, and one approximate scheme, the 
Fixman-Freire approximation, for computing DNA melt¬ 
ing profiles. 

Let us first introduce some notation. The full sequence 
of basepairs enters via the position-dependence of the 
statistical weights 

ai = exp{-F;hb(i)/[i?T]} (Al) 

for breaking the hydrogen-bonds of the basepair at posi¬ 
tion i, and 

S, = eM-Estii-l,i)/[RT]} (A2) 

for disrupting the stacking interactions between basepairs 
i — 1 and i. In addition, we need the parameters ^ and 
the loop factor g{m), see main text. 


1. The Poland algorithm 

We now provide an exact scheme, the Poland 
algorithnP^ for clamped ends, for computing the proba¬ 
bility that basepair k is open. The algorithm presented 
here extends previous methods by incorporating explicit 
stacking interactions: previous approaches “lumped” one 
Boltzmann factor of stacking interactions together with 
the giving rise to a cooperativity parameters, (Tq. 

The Poland algorithm has two components. First, the 
conditional probability Pc{k) that basepair A:-I-1 is closed, 
provided that basepair k is closed, is computed. Note 
that we here use the (somewhat funny) short-hand nota¬ 
tion Pc(fc) = Pc{k+l\k) following Poland.!^ Secondly, the 
unconditional probability Pu{k) that basepair k is closed, 
is computed. In practice, however, one does not directly 
work with Pdk) but rather introduces = UkSkPcik). 
If we now extend the algorithm from Ref. H^to the case 
of clamped boundary conditions and explicit stacking pa¬ 
rameters, we find the following recursion relation: 

Tn = 

N-k 

n = akS,d+J2 k = N-l,N-2,...,0, 

i=i 

(A3) 

where a/j(j) = qi{k, j)Yli^k+i EM) with the statistical 
weight qi{k,j) for a bubble with its left boundary at 


position k (basepair k is closed and basepair fe -|- I is 
open) and consisting of j consecutive open basepairs: 

j+k j+k+l 

qi{k,j) = MU) nn <5/. (A4) 

/ —fc + 1 Z —fc + 1 


We above defined Oq = i5o = 1. The quantity atU) is 
conveniently computed using the recursion relation 

afe(l) = MU)^k+2Tk+i k = 0,l, N -1, 
akU) = f.^^\M k+iak+iU - 1) j = 2,...,N -k. 

gU - 1 ) 

(A5) 

The (unconditional) probability P«(fc) that basepair k is 
then obtained from the equations: 


Pu{0) = I, 

PuU)=To, 

k-l 

Pu{k + 1) = Pu{k)Pdk) + Y, Mj)PuU), fc = 1,2,..., N. 

j=o 

(A6) 

We further define fik{j) = QiUM — j)]^.^ PM) from 
which we find the recursion relation: 

Mk - 1) = MU) Tk-iTk /c = 2, 3,..., N, 

Oik-lOk-l 

/.X 9{k-j) Sk+lrj. 

i = 0,I,...,fc-2. (A7) 


The Poland algorithm for clamped DNA and explicit 
stack ing parameters thus becomes: First use Eqs. (A3) 
and (A5) to obtain Tj. for k = N,N — 1,...,0. Then 
use Eqs. (A6) and (A7| to obtain Pu{k) for k = 

0,1,., N + 1. A consistency check of the calculation 

is that P„(V -I- I) = 1 due to the clamping. 


The mean number of bubbles used in appendix |G l| can 
easily be determined from P„(fc) and P^k) as the number 
of left boundaries (basepair i closed, basepair i + 1 open): 


Mb) = 

i 


= Y “ XI PM)PM) 

i i 


(AS) 


where we introduced the “spin” of a basepair: Si = I if 
basepair i is closed and Si = 0 otherwise. 

The Poland algorithm scales with the number of base- 
pairs N as In practice, we therefore do not use the 
Poland algorithm for V > 10 kbp. In the next section we 
present an algorithm whic scales linear in V, the Fixman- 
Freire approximation.!^ 
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2. Fixman-Preire approximation 


The Fixman-Freire (FF) approximation uses the fact 
that a sum of exponentials is a good approximation for a 
power law function. The sums in Eqs. (A31 and (A6) give 
rise to the N'^ scaling of the Poland algorithm. Replacing 
the loop factor g{m) (as part of both sums, see Eq. (A4)) 
by a sum of exponentials with I terms, reduces the scaling 
of the approximation to an A^I-scaling. 

Thus, within Fixman-Freire’s method, the following 
approximation is introduced: 


basepairs, basepairs number 0 and A^ + 1, are here as¬ 
sumed clamped (closed). Garel and Orland, see Ref. [HI 
proceed by deriving a recursion relation (see below) for 
the statistical weights (partition functions), and 

(5fl(k), for the “left” and “right” regions with the con¬ 
straint that basepair k is closed. Given QlO^) and Qn^k) 
we can write most quantities of interest. In particular, 
the probability for having basepairs k closed is simply 


P{k) = 1 - 


QLik)QR{k) 

Ql(a^+1) ■ 


(Bl) 


gim) ^'^Ci exp{-bim) 


(A9) 


. Applied on the sum in Eq. (A3), we find 


Tk = akSk[l+ '^CiSi{k)] ^ 


(AlO) 


and the recursion relation 
S,{N) = 0 

Si{k - 1) = Tfeexp(-6i)[^4+i + -S'j(A:)]. (All) 


Similarly, Eq. (A6) becomes 


Due to symmetry, recursion relations for QR{k) are of 
the same form as those for QL{k). We therefore now 
limit ourselves to studying QL(k). We further find it 
convenient to rescale (5i,(fc) as 


.. N _ Qdk) 

<lL{k) = 


Pif32 • • • /3fc ’ 


(B2) 


where Pk = 0:^4 with and defined in Eqs. (Al) 


and (A2). In terms of these rescaled statistical weights, 
the Garel-Orland recursion relatiorfi^ becomes: 


fc-i 

Pk+iqL{k + 1) = qL{.k) + cro,fe+i ^ g{k - j)qL{j) (B3) 

1=0 


Pu(,k + 1) = Puik)P,ik) + ^ c,A{k) (A12) 

i=l 

with recursion relation 

A,(0) = 0 

Ai{k -I-1) = Pc{k + l)afc+i4+i exp(-6i) 

[^6k+iPu{k)Pcik) + Mk)]. (A13) 

The coefficients Ci and bi are determined from a fit on 
21 points, li (chosen such that the log’s are equally dis¬ 
tributed on the domain where the approximation should 
be valid, i.e. from 1 to the total length of the DNA se¬ 
quences). Starting with all coefficients set to zero, first 
the guess values c/ and bj are determined from gihi) 
and g(l 2 i-i), then follow ci-i and 6/_i using the new 
set of parameters and so on. If no convergence can be 
achieved, the parameters are determined for a larger in¬ 
terval, which is then slowly shrunk again with the new 
set of parameters as a first guess. 


Appendix B: DNA melting behavior for homo-DNA 


for k = 1, We above introduced the local coopera- 

tivity parameter cro^fe = ^4- The “boundary” conditions 
are qL{0) = 1 and (7 l( 1) = 1/4- We note that, in prin¬ 
ciple, the recursion relation above can be used to predict 
DNA melting curves. However, the quantity qL^k) typi¬ 
cally “explodes” exponentially with the number of base- 
pairs. Therefore, two alternative numerical methods, the 
Poland algorithm and the Eixman-Freire approximation, 
are described in appendix [A| 

Let us consider the case of homo-DNA, i.e. when the 
basepair energies are the same along the DNA. We then 
set Pk = P and tro.fe =_£o in Eq. (B3). We notice that 
the last term in Eq. (B3) has the form of a convolution 


(if we add the term j = k to the sum). We therefore 
introduce the generating function (if we replace z —>■ l/z 
below we get a z-transform) 


9L 


iz) = ^qUk)z^ 


(B4) 


k=0 


which, after using the convolution theorem for z- 
transforms and some rearrangements (for convenience we 
also define g(0) = 0), gives the solution to Eq. (B3| as: 


In this appendix we analyze exact recursion relations 
for the partition function for the Poland-Scheraga model, 
and provide an analytic expression for the melting prob¬ 
ability for homo-DNA. 

Eollowing Ref. [T8| we divide the DNA molecule into 
a left region, i.e. basepairs {1,..., k — 1} and a right re¬ 
gion, basepairs {k + 1,..., A}, where k varies. The end 


Ql{z) = - -(B5) 

1-/3 ^z[l + aog{z)] 

where 

OO 

9i.z) = '^g{k)z^ (B6) 
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FIG. 9. Contours used for solving the Poland-Scheraga model 
analytically for homo-DNA. Integration along the black cir¬ 
cle around 2 = 0 gives the statistical weight qL{k) through 
Eqs. (B5l and (B7l. The blue contours depict a useful de¬ 
formation of the black contour: We deform the contour to 
“infinity”. When deforming the original contour, we “pass” 
the pole at 2 = 2o, which must be subtracted (notice the 
“negative” orientation of the contour around 2 = 2o). The 
contour along Re 2 = [1, oo] is needed for the case of a power 
law form for the loop factor, as the associated 2 -transform 
(polylogarithmic function) has a branch cut there. 


For the case g{k) = we have that g{z) = Lic( 2 ) 
is the polylogarithmic function. It is noteworthy that 
analytic continuation of the polylogarithmic function to 
the complex plane has a branch cut for Re 2 > 1. The 
generating function qL{z) therefore also has a branch cut 
in the same domain. 


In order to invert the exact result for the generating 
function, Eq. (B5), we note, see Eq. (B4|, that qL{k) is 
the /c:th term in a series expansion in 2 , i.e., 


qhik) 


1 d’^qdz) . _ ^ ^ 

fc! dz^ 27rz 


<1l{z) 


dz 


(B7) 


where the last equality is a standard result in complex 
analysis.l^ The integration above is over a closed con¬ 
tour C around 2 = 0, se e Fig . Following Ref. ITCIthe 
contour integral in Eq. ( |B7[ ) is evaluated by deforming 
the original contour as described in Fig. The deformed 
contour has several parts. One of these parts is a contour 
around the pole zq of qL{z), and, in fact, the thermody¬ 
namic behavior is determined by the contribution from 
this pole.^ For large k we therefore have 


(B8) 

where Cq is the contour around the pole depicted in Fig. 
The pole, 20 , is determined by the condition 9 ^( 20 ) = 


0 , i.e.. 


9{zo) = 


Id- zq 


(B9) 


where we used Eq. (B5). So, using the Cauchy theorenP^ 
we have: 

^-k-l 


qL{k) ~ -Res[- 


,2 = 20 ] (BIO) 


1-/3 i2[l -b aog{z)]'' 

where Res[...] denotes the residue. Using the fact thalP^ 
Res[A( 2 )/i 3 ( 2 ), 2 = 2 o] = A{zo)/B'{zo) if zq is a simple 
pole and A{z) is a regular function, we find: 


dLik) = 




1 


(Bll) 


CTo[s'(2o) + zog'{zo)]' 

This result holds provided that 20 G [0,1]; if 20 > 1 then 
the pole is in a branch cut (for a power-law form of g{m)), 
see Fig. and does not contribute toq^lk). Finally, we 
simplify t he re sults above using Eq. (B9), and combining 
with Eq. 


we find the partition function: 

— k 

^ 1 -b CTo2o/3-1/i(2o) 


(B12) 


where we defined the function: h(z) = zg'(z). Eq. (B12) 


is our final expression for the rescaled partition function 
for the Poland-Scheraga model for homo-DNA melting 
in the thermodynamic limit. For the general case we 
evaluate g(z) using the definition, Eq. (B 6 I. Using this 
same definition, we find that 


= X! ^9{k)z'^ 


(B13) 


k=l 


thereby providing an explicit expression for all quantities 


in Eq. (B12|). 
Final 


y, by combining Eqs. (Bl), (B2), and (B12), we 


obtain our final expression for the melting probability: 

1 


P = 1 - 


l-bcro/3 ^ 20 / 1 ( 20 ) 


(B14) 


where g[z) is given in Eq. (B 6 I and h{z) in Eq. (B13|. 


These expressions can be used to numerically evaluate 
g{z) and h(z) by replacing the series by sums. The 
quan tity Zq is determined by (numerically) solving Eq. 
(B9). Thus, the equations given here provide means 


for straightforward calculation of the melting probabil¬ 
ity for different choices of g{m) in the thermodynamic 
limit. For the case g{k) = k~'^ we have g{z) = Lic( 2 ) 
and h{z) = Lic_i( 2 ), where Lic( 2 ) is the polylogarithmic 
function. 


Appendix C: Analysis of melting curve for 
unconfined DNA 

This appendix provides results for melting curves of 
unconfined DNA, and, thus, serves a background to the 
main text, where melting of confined DNA molecules is 
studied. 
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1. Deviation from the Marmur-Doty melting 
temperature 


The Marmur-Doty formula, Eq. , is based on a 
simple estimate where one assumes melting to occur at 
a temperature where the free energy for breaking one 
hydrogen bond and one stacking interaction equals zero. 
Due to, for instance, neglect of the “boundary energy” 
(l/2)iiTlog^ and the loop factor g{m), which favors 
smaller bubbles, Eq. Q is only an approximation. 

In order to estimate shifts in melting temperature com¬ 
pared to the Marmur-Doty formula for unconfined DNA, 
we note that an entropic contribution i?log^ has to be 
assigned to each bubble. The corrected melting temper¬ 
ature can then be estimated from the number of bubbles, 
Nb, as 


rr^B _ 

— 


EH 

Es 


'ref 


Sre!-mOg^^ 


Tm 


(Cl) 


In practice, we further replace Nb by its expected num¬ 
ber of bubbles {Nb)- In Eigure[l0]we tested the above 
mean-held prediction in the following way: using the 
Eixman-Ereire approximation we determined the shift in 
melting temperature from the Marmur-Doty formula for 
a range of pat and loop exponents c. Simultaneously, 
we determined the expected number of bubbles, (IVs), 
see Sectio n [A] f or details. From Figure [T0| we notice that, 
while Eq. (Cl) captures the typical behavior for the shift 
as a function of number of bubbles, there are deviations. 
For c ^ 1, there is clear relation between the number of 
bubbles and the deviation from the Marmur-Doty melt¬ 
ing temperature. For larger values of c, however, the 
trend is less pronounced. A possible explanation for this 
“large c” effect is that while the curves get wider with 
smaller values of c, the asymmetry of the melting tran¬ 
sition results in the point T{f = 0.5) to be shifted to 
smaller temperatures. 

Notice that we in the main text avoid “problems” 
with the Marmur-Doty formula by always using simu¬ 
lated melting temperatures as reference. 

A final word of caution is here in order. The hydro¬ 
gen bond energies we use in the simulations are not mea¬ 
sured for individual base-pairs, but estimated from global 
melting of a few model sequences and therefore already 
contain the contribution of the ring factor. We therefore 
get too large temperatures from our simulations already 
for unconfined DNA with c = 2.12, wh ere th e Marmur- 
Doty relation was found experimentally^i^Ell For smaller 
values of c (or smaller channels), the melting transition 
gets broader and the mean number of bubbles at Tm in¬ 
creases. 


2. Width of the melting transition as a function of 
loop exponent c 

Figure [TT] shows the change of the width a for different 
loop exponents c for unconfined DNA. The width was ex- 
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FIG. 10. Shift of the simulated melting temperature rM,sim 
(units of Kelvin) compared to the Marmur-Doty melting tem¬ 
perature Tm.m.d. as a function of the simulated mean num¬ 
ber of boundaries at TM.aim at TM.aim for random uncon¬ 
fined DNA. The solid line shows our theoretical estimate 
from Eq. (Cl I. For c ^ 1, there is clear relation be¬ 
tween the number of bubbles and the deviation from the 
Marmur-Doty melting temperature. For larger values of c, 
however, we find less correlation between the shift in melt¬ 
ing temperature and {Nb}- Averages are taken over 200 
different sequences of length 2 • 10® basepairs at AT ratios 
Pat = 0,0.25,0.5,0.75,1 (different markers), each for differ¬ 
ent loop exponents c = 2.12,1, 76,1.5,0.5, 0 (indicated by dif¬ 
ferent colors). 


tracted using the same approach as described in subsec¬ 
tion |TV^ We find that increasing c decreases the width 
of the transition as it should. These findings hold for the 
three cases of interest here: 1. homo-DNA, 2. random 
DNA for different AT fractions, and 3. T4 phage DNA. 
The prefactors for the width are larger for case 2. than 
for case 1. The homo-DNA case does not give identical 
results to the /at = 0 and /at = 1 cases due to the ad¬ 
ditional different stacking parameters for the latter cases 
as compared to the homo-DNA case, see section m 


3. Analytical prediction for different loop 
exponents c 


Figure [T^ shows the melting probability for unconfined 
homo-DNA for different values of the loop exponent c, 
both using the numerical Poland algorithm and the an¬ 
alytic prediction introduced in section As before, we 
used Eq. ( B14 ), where Zq was numerically obtained by 
solving Eq. (B9). We find very good agreement between 
the numerical and analytical result. Note, however, that 
due to the finite size of the DNA sequences the numeri¬ 
cal result for c = 1.5 deviates slightly from the analytical 
result, which is obtained for the thermodynamic limit. 
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♦ numerical 
T analytic estimate. 
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FIG. 11. Width of the melting curve as a function of the loop 
exponent c for unconfined homo-DNA (top), unconfined ran¬ 
dom DNA (middle) and unconfined T4-phage DNA (bottom) 
(with associated random DNA, same AT fraction). Parame¬ 
ters are the same as in Fig. 



-c= 1.5 
-c = 0.5 
c = 0 
analytic 


FIG. 12. Melting probability as a function of temperature for 
an unconfined homo DNA for different loop exponents. The 
numerical results for a sequence of length 50 kilo basepairs 
(solid lines) are compared to the analytic prediction in the 
thermodynamic limit (red marks). We used ^ = 10® and 
melting temperature for a sequence ‘A ... A’, see Eq. (§. 
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